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Abstract 

At high-energies the gluon-gluon fusion is the dominant mechanism of cc production. This 
process was calculated in the NLO collinear as well as in the k^-factorization approaches in the 
past. We show that the present knowledge of gluon distributions does not allow to make a precise 
predictions for cc production at LHC, in particular at forward rapidities. In this paper we study 
production of cc pairs including several subleading mechanisms. This includes: gg — )• QQ, jg — >■ 
QQi 91 ~^ QQi 77 ~^ QQ- III this context we use MRST-QED parton distributions which 
include photon as a parton in the proton as well as elastic photon distributions calculated in the 
equivalent photon approximation. We present distributions in the c quark (c antiquark) rapidity 
and transverse momenta and compare them to the dominant gluon-gluon fusion contribution. We 
discuss also inclusive single and central diffractive processes using diffractive parton distribution 
found from the analysis of HERA diffractive data. As in the previous case we present distribution in 
c (c) rapidity and transverse momentum. Finally we present results for exclusive central diffractive 
mechanism discussed recently in the literature. We show corresponding differential distributions 
and compare them with corresponding distributions for single and central diffractive components. 
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I. INTRODUCTION 



In the past we have calculated inclusive cross section for heavy quarks production at 
hadron colliders. These calculations wereperformed using an approach based on the un- 
integrated parton distributions functions [l|, |2|- It is known that gluon-gluon fusion is the 
dominant mechanism at high energy. However, other mechansims were not carefully studied 
in the literature. 

It is the aim of this work to present contributions of several subleading terms usually 
neglected in the analysis of cc production. We wish to include contributions of photon-gluon 
(gluon-photon) as well as purely electromagnetic contributions of photon-photon fusion. 

We wish to discuss also diffractive processes (single and central) in the framework of 
Ingelman-Schlein model corrected for absorption. Such a model was used in estimation of 
several diffractive processes 

The absorption corrections turned out to be necessary to understand a huge Regge- 
factorization breaking observed in single and central production at Tevatron. 

Recently a surprisingly large cross section for exclusive cc production has been reported 
joj. Her we will show results for RHIC and LHC energies. 



II. PRODUCTION OF HEAVY QUARKS 

In the leading-order (LO) approximation within the coUinear approach the quadruply 
differential cross section in the rapidity of Q (yi), in the rapidity of Q (2/2) and the transverse 
momentum of one of them [pt) can be written as 

dcr ] 



dyidy2d?Pt IGvr^s^ 



^a;iPi(xi,/x^) X2Pj{x2,n^) \M,j^qq\^ . (2.1) 



Above, pi(xi,/i^) and Pj(x2,/i^) are the familiar (integrated) parton distributions in hadron 
hi and /i2, respectively. There are two types of the LO 2 — )■ 2 subprocesses which enter 
Eq.f l2.ip : gg — > QQ and qq — )■ QQ. The first mechanism dominates at large energies and 
the second one near the threshold. In particular for the gluon-gluon fusion the cross section 
formula takes a simple form: 

da 



There are three (s, t and u) diagrams in the leading order 10 | 



The parton distributions are evaluated at: xi = ^ (exp(|/i) -|- exp(?/2)), X2 



^ (exp(— ?/i) + exp(— ^2)); where rrit = ^pi + mq. The formulae for matrix element squared 
averaged over the initial and summed over the final spin polarizations can be found e.g. in 
Ref.^OI]. 

The inclusive heavy quark/ ant iquark production can be also calculated in the framework 
of the fct-factorization. In this approach transverse momenta of initial partons are included 
and emission of gluons is encoded in so-called unintegrated gluon distributions (UGDFs) 

In the leading-order (LO) approximation within the fcrfactorization approach the quadru- 
ply differential cross section in the rapidity of Q (yi), in the rapidity of Q (1/2) and the 
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FIG. 1: Standard diagrams representing mechanisms for production of heavy quarks. 



transverse momentum of Q {pi,t) and Q {p2,t) can be written as 

/ — — —T7^Z27:7zr:A^ \^^j^QQr (2-3) 



dyidy2d'^Pi^td'^P2,t 7^ ^ ^ lQn'^{xiX2sy 

6"^ + K2,t - Pi,t - P2,t) J^i{xi, J J^j{x2, ni t) , 

where J-i(xi,K^J and J-'j(x2, k^.J are so-called unintegrated gluon (parton) distributions. 
Now the unintegrated parton distributions must be evaluated at: xi = ^^exp(?/i) + 

^exp(?/2), X2 = ^exp(-?/i) + ^exp(-?/2), where = ^plt + "^q- 

III. PHOTON INDUCED PRODUCTION OF HEAVY QUARKS 
A. MRST-QED parton distributions 

As discussed above the dominant contributions are initiated by gluons or quarks and 
antiquarks. In general even photon can be a constituent of the proton. This was considered 



only in one work by Martin, Roberts, Stirling and Thorne [Ul]. Below we repeat the main 
aspects of their formalism. 

The factorization of the QED-induced collinear divergences leads to QED-corrected evo- 



lution equations for the parton distributions of the proton [11 

dqi{x,n'^) as dy r X 2^ , -p f ^ 2^\ 

2n Jx y ^ y y ' 

dg[x,n'^) as dy r , ,x 2n , p / ^\\ 

9.(x,^=) _ -j^i^lp^^^^) Y.^),,(-,^f) + P-Ml(-,^')}, (3.1) 

zn Jx y ^ ■ y 17 



d log /i^ 2tc Jx y 



where 



P — r*"! P P — r^^ P 
nn ^ F <i<ii iQ F gg^ 

1 2 2 
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FIG. 2: Diagrams representing mechanisms for production of heavy quarks, which included pho- 
tons. 



and the parton distributions fullfil momentum conservation: 

(ia: X I + 5f(x,/i^) + 7(x,/i^)| = 1 . (3.2) 



B. Mechanisms w^ith one or two photons 

If the photon is a constituent of the nucleon then other mechanisms presented in Figl2] 
are possible. 

Here the cross section can be calculated similarly as for the gluon-gluon fusion. A corre- 
sponding triple differential cross section can be written as: 

1 



a;i7m(a;i,/x ) X2g{x2,n ) {M^g^Qgl 



dyidy2d?Pt IGvr^s^ 



Xi7i„(Xi,/i ) X21in{x2,^l ) IA^^^^qqP (3.3) 



for gluon-photon, photon-gluon and photon-photon contributions, respectively. Compared 
to gluon-gluon case here only t and u diagrams occur. 

The above contributions include only cases when nucleons do not survive a collision and 
nucleon debris is produced instead. The case when nucleon survives a collision has to be 
considered separately. In this case one can include corresponding photon distributions where 
extra "el" index will be added to denote that situation. Corresponding contributions can be 
then written as: 

c/a^T-' 1 



dy,dy2d-p, levr^J^-^^^^^^' ^ ^ ^ ^ l-^- 

da^-'9 1 



dy,dy2d-p, - 16^^^^^'^^^'^ ) ^^^(^^'^ ) l-^7.->QQr , 
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dyidy2(Ppt 




da"''''"'''' 
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dyidy2d?Pt 





XlJin{Xi, fi"^) X2jel{x2, fX^) \M 



77-s>QQl 5 



Xi-fel{Xl, jJ, ) X2%n{x2,IJ' ) 



77-s>QQl 5 



Xi7ei(Xi,/i ) X2-iel{x2,fJi ) 



77— s>C 



(3.4) 
(3.5) 



The elastic contributions are calculated using Drees-Zepenfeld (elastic) parametrizations of 
photon fluxes 12| which include nucleon electromagnetic form factors. 



IV. RESULTS 



A. Gluon distributions and small-x region and its relation to cc production 



In Figl3]we show three different leading-order gluon distributions from the literature 11 



13l . Il4l | (left panel) and photon distributions ll| (right panel) as a function of longitudinal 
momentum fraction x for a fixed scale /i^ = 10 GeV relevant for cc production. Above x > 
10~2 all the distributions coincide. For smaller values of x they diverge and can be different 
by almost an order of magnitude. What are consequences of this divergence for cc pair 
production? This will be discussed below. 
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FIG. 3: Different leading order gluon distributions form the literature for the factorization scale: 
/i^ = 10 GeV^ (left panel) and leading order photon distributions for factorization scale: /z^ = 10 
GeV^ (right panel). 

Before we go to cross sections in transverse momentum and rapidities, in FigJUwe present 
distribution of the cross section in ^1 = logioXi and ^2 = logioa;2 for two different energies 
a/s = 500 GeV (updated RHIC) and a/s = 14 TeV (nominal LHC energy). One can clearly 
see that the Xi and X2 values are strongly correlated. Typical values at ^/s = 500 GeV are 
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xi,X2 ~ 0.5 • 10~^ and at ^/s = 14 TeV are Xi,X2 ~ 10"'*. In the latter case x's as small as 
10~^ may appear in the forward c or c region. This is clearly a region of x which was never 
studied so far. 




FIG. 4: Distributions in xi and X2 for two different energies: ^/s = 500 GeV (left) and ^/s = 14000 
GeV (right). In this calculation GRV94 gluon distributions have been used. 



Now let us present distributions in transverse momentum of c (or c) for gluon-gluon 
fusion mechanism for different gluon distributions and different popular choices of scales 
(yU^ = 4ml, invariant mass of the cc system M^g, +m^). We show our results for ^/s = 500 
GeV (Figl5]) and a/s = 14 TeV (Fig|6]). One can clearly see that for some choices of gluon 
distribution function and scales the results for ^/s = 14 TeV are not physical. This shows 
how badly known are gluon distributions at the low x. 




FIG. 5: Distribution in quark/antiquark transverse momentum at ^/s = 500 GeV for different 
choices of scales and for different gluon distributions: GRV94 (left panel), MRST2004 (midle 
panel) and MSTW2008 (right panel). In this calculation we have used fxj^ = = ^roq. 
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FIG. 6: Distribution in quark/antiquark transverse momentum at ^/s = 14 TeV for different choices 
of scales and for different gluon distributions: GRV94 (left panel), MRST2004 (midle panel) and 
MSTW2008 (right panel). In this calculation we use have used = /u^ = 4mQ. 



B. 'jg and subprocesses 

In FigJT] and in FigJH we show results for different gluon distribution functions for the 
RHIC energy ^/s = 500 GeV and nominal LHC energy ^/s = 14 TeV, respectively. At the 
LHC energy the results for different GDFs differ considerably which is a consequence of 
the small-x region as discussed in the previous section. The differences at the nominal LHC 
energy ^/s = 14 TeV are particularly large which can be explained by the fact that a product 
of gluon distributions (both at small x) enters the cross section formula. A new measurement 
of cc at the nominal LHC energy will be therefore a severe test of gluon distributions at 
small X and not too high factorization scales not tested so far. Similar uncertainties for the 
jg and (77 are smaller as here only one gluon distribution appears in the corresponding cross 
section formula. 

It is very difficult to quantify uncertainties related to the photon PDFs as only one set of 
PDFs includes photon as a parton of the proton. Here the isospin symmetry violation (not 
well known at present) would be an useful limitation. Our collection of the results for the 
photon induced mechanisms show that they are rather small and their identification would 
be rather difficult as the different distributions are very similar to those for the gluon-gluon 
fusion. Our intension here is to document all the subleading terms in one publication. Our 
etimation shows that the sum of all the photon induced terms is less than 0.5 % and is by 
almost 2 orders of magnitude smaller than the uncertainties of the dominant leading-order 
gluon-gluon term. 
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FIG. 7: Transverse momentum distribution for the standard gluon-gluon mixed gluon-photon and 
photon-gluon as well as for photon-photon contributions for RHIC (left panel) and LHC (right 
panel). Three different gluon distributions were used. The photon distributions are from We 
show contributions when proton survives the collision (called elastic) and when hadronic debris is 
produced (called inelastic). 



V. SINGLE AND CENTRAL DIFFRACTION 



A. Formalism 



The mechanisms of the ordinary as well as diffractive production of heavy quarks (cc) 
are shown in Figs l9|T0l 

In the following we apply the Ingelman and Schlein approach^. In this approach one 
assumes that the Pomeron has a well defined partonic structure, and that the hard pro- 
cess takes place in a Pomeron-proton or proton-Pomeron (single diffraction) or Pomeron- 
Pomeron (central diffraction) processes. We calculate triple differential distributions as 



00 



dyidy2dpi 



K 



M 



(5.1) 



^ In the literature also dipole model was used to estimate diffractive cc production 15 1 
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FIG. 8: Different contributions to distributions in rapidity of c quark/antiquark at ^/s = 14 TeV 
for different gluon distributions: GRV94 (left panel), MRST2004 (midle panel) and MSTW2008 
(right panel). In this calculation we have used fip = iJ-j^ = s. 





FIG. 9: The mechanism of single-diffractive production of cc. 



da 



SD 



K- 



M 



dyidy2dpt IGvr^s^ 



(5.2) 



da, 



CD 



K- 



M 



dyidy2dpt IGn'^P 



(xiqf (xi, /i^) X2qf {x2, /i^)) + (xiqf (xi, /i^) X2qf {x2, 



(5.3) 
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FIG. 10: The mechanism of central-diffractive production of dileptons. 



for ordinary, single-diffractive and central-diffractive production, respectively. 

We do not calculate the higher-order contributions and include them effectively with the 
help of a so-called iiT-factor. We have checked that this procedure is precise enough in the 
case of ordinary Drell-Yan process. The i^'-factor is calculated as for the Drell-Yan process 



Here the running coupling constant = ^^(/i^) is evaluated at /i^ = M^q. 

The 'diffractive' quark distribution of flavour / can be obtained by a convolution of the 
flux of Pomerons /ip(xip) and the parton distribution in the Pomeron gj/]p(/3, /i^): 



/p dxjp Oil 
dxjpdf36{x-xjp(3)qf/-sp{/3,f/)fjp{xjp) = / /ip(a;ip)g//ip( — 
Jx Xjp Xjp 

(5.4) 

The flux of Pomerons /ip(a;ip) enters in the form integrated over four-momentum transfer 

f ^max 

fjp{xjp) = I dtf{x-sp,t), (5.5) 



with tmin,tmax being kinematic boundaries. 

Both pomeron flux factors f-sp{x-p,t) as well as quark/antiquark distributions in the 
pomeron were taken from the HI collaboration analysis of diffractive structure function and 
diffractive dijets at HERA 1^. The factorization scale for diffractive parton distributions 
is taken as /i^ = s. 



B. Results 

Let us start presentation of our results for diffractive mechanisms. 

In Fig{TT]we show transverse momentum distributions of charm quarks (or antiquarks). 
The distribution for single diffractive component is smaller than that for the inclusive gluon- 
gluon fusion by almost two orders of magnitude. Our results include gap survival factor. 
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Corresponding values are taken the same as in Ref. 17|]. The cross section for inclusive 
central diffractive component is smaller by additional order of magnitude. In addition we 
show the cross section for fuUy^ exclusive mechanism discussed in section IVII Below we shall 
use the following notation: 00 for standard nondiffractive component, Od or dO for single 
diffractive and dd for central diffractive components. 
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FIG. 11: Transverse momentum distribution of c quarks (antiquarks) for RHIC energy ^/s = 500 
GeV for three different parton distributions. The result for single diffractive (Od or dO), central 
diffractive (dd) mechanisms are compared with the standard gluon-gluon fusion contribution (00). 




FIG. 12: Rapidity distribution of c quarks (antiquarks) for RHIC energy ^/s = 500 GeV for three 
different parton distributions. The result for single diffractive (Od or dO), central diffractive (dd) 
mechanisms are compared with the standard gluon-gluon fusion contribution (00). 



^ Although the calculation assumes simple cc state hadronization leads to more complicated states 
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In FigJTS] we show similar results for nominal LHC energy ^/s = 14 TeV. The situation 
and the interrelations between different components is qualitatively the same. Here some- 
what smaller gap survival factors were used 17|. The distributions for all components are 
somewhat broader than those for the RHIC energy shown above. 
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FIG. 13: Transverse momentum distribution of c quarks (antiquarks) for RHIC energy ^/s = 14000 
GeV for three different parton distributions The result for single diffractive (Od or dO), central 
diffractive (dd) mechanisms are compared with the standard gluon-gluon fusion contribution (00). 



In FiglHlwe show distributions in quark (antiquark) rapidity. We show separately contri- 
butions of two different single- diffractive components, which give the same distributions in 
transverse momentum in FigJT3l When added together they produce a distribution similar 
in shape to the standard inclusive case. Here different parton distributions functions give 
similar result. The distributions for different proton gluon distributions are quite different. 
This was already observed when discussing photon induced components in section IVI[ 

Also two dimensional distributions can be interesting as here different mechanisms may 
occupy different parts of the phase space. 

In Fig{T5]we show distributions in the rapidity of the pair and quark-antiquark invariant 
mass. Although the distributions are somewhat different the differences occur in the regions 
which may be difficult to measure. The spread in the pair rapidity for the central diffractive 
component is much smaller then that for the inclusive case. 

Finally we show distributions in quark and antiquark rapidities. The distribution for the 
inclusive central diffractive mechanism are concentrated at midrapidities. This is a rather 
universal feature of diffractive processes. 

The cross section for single and central diffraction is rather small. However, a very 
specific final state should allow for its identification by imposing special conditions on the 
one-side (single-diffractive process) or on both-side (central diffractive process) rapidity gaps. 
We hope that such an analysis is possible at LHC. Special care should be devoted to the 
observation of the exclusive cc production where the observation of D mesons associated 
by a few pious would be a proper signal [isj. Without a special analysis of the final state 
multiplicity the exclusive cc production may look like an inclusive central diffraction. At 
present there is no analysis of the final state production for the exclusive cc. 

A comparison of cross sections for both components will be done in the next section. 
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FIG. 14: Rapidity distribution of c quarks (antiquarks) for LHC energy y/s = 14 TeV for three 
different parton distributions. The result for single diffractive (Od or dO), central diffractive (dd) 
mechanisms are compared with the standard gluon-gluon fusion contribution (00). 



VI. EXCLUSIVE CENTRAL DIFFRACTIVE PRODUCTION OF cc 

There is recently a growing theoretical interest in studying central exclusive mechanisms 
of different particles production at high energies, which constitute a special category of 
double-diffractive processes. To date, only a few exclusive processes have been measured so 
far at the Tevatron collider (see [l^ and references therein). In particular, central exclusive 
production (CEP) of the Higgs boson is a flag process of special interest and importance in 



the upcoming Higgs searches at the LHC (see e.g. Ref. |2d. l21|). 

Generally, in the case of the central exclusive production pp — )■ pXp with the leading 
protons, the central system X should necessarily be produced in the color singlet state, such 
that the proton remnants and the X system are disconnected in the color space and their 
hadronisation occurs independently giving rise to rapidity gaps [22| . From the experimental 
point of view, CEP procesess are very attractive, because of the rare clean experimental 
environment, related to = selection rule, and great mass resolution of the centraly 
produced object. Such unique features give a new possibility to exploit hh high branching 
ratio decay channel of the Higgs boson, which is rather impossible in standard inclusive 
measurements, due to very large QCD background. Therefore, the QCD mechanism of 
central exclusive heavy quark dijets is a source of the irreducible background to the exclusive 
Higgs boson production. 

Central exclusive production of cc and bb pairs was studied in detail in our previous 
papers jo], 2^ 21|. In these calculations the pp — )• p{qq)p reaction, illustrated in Fig. [T71 was 



considered as a genuine 4-body process with exact kinematics. The applied perturbative 
model of theoretical predictions is based on the Khoze-Martin-Ryskin (KMR) approach 
used previously for the exclusive Higgs boson production [Isj. Total cross sections and 
differential distributions for heavy quarks are calculated by using fct-factorization approach 
with help of the KMR unintegrated gluon distribution functions. 

This QCD model works very good in the case of exclusive dijets and charmonia produc- 



tion, what was confirmed by CDF data |24j-|27l|. However, estimated uncertainties related 
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FIG. 15: Two-dimensional distributions in rapidity of the pair and the quark-antiquark invariant 
mass for standard (upper left), single diffractive (upper right and lower left) and central diffractive 
contributions. In this calculation MRST04 distributions were used. 



to gluon densities, factorization and renormalization scales, as well as due to absorption 
corrections are quite large. It makes the situation somewhat clouded and prevents definite 
conlcusions, especially in the case of the exclusive production of heavy quark pairs. In this 
context, the promising idea, how to clarify and calibrate purely- known parameters of the 
theoretical model, is to study cc cross section by exclusive measurements of DD meson pairs. 
Such experimental studies are being performed now at Tevatron and could be also available 
in Run II experiments at RHIC and at LHC 

Therefore, it is also very interesting, from both, theoretical and experimental side, to 
compare mechanism of central exclusive production of charm quarks with standard single 
and double diffractive processes. Such an analysis of differential cross sections has never 
been done before but could bring important informations about differences in kinematics 
and in production rates between them, what is crucial for future measurements. 



According to the KMR approach [23|, |28|, |29j we write the amplitude of the exclusive 
diffractive qq pair production pp — )■ p{qq)p as 

s 7r%e2 ^ r ff^{xux[, ql^, qf ^, ti)/g°f (xa, 4, qlt, qjt, ta) 



-^A,A, = ■ 5> / d\o,t V^lll{q„ q,) X 



qlt qlt qlt 
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FIG. 16: Two dimensional distribution in rapidity of the quark and rapidity of the antiquark 
for standard (upper left), single diffractive (upper right and lower left) and central diffractive 
contributions. In this calculation MRST04 distributions were used. 



Pi pi 




FIG. 17: The mechanism (lef panel) and kinematics (right panel) of exclusive double-diffr active 
production of heavy quarks. 



where Xg, Xg are helicities of heavy q and g, respectively, ti_2 are the momentum transfers 
along each proton line, qi^t, (l2,t, a^i,2 and qo^t, x[ ~ x'a ^ Xi^2 are the transverse momenta and 
the longitudinal momentum fractions for active and screening gluons, respectively. Above 
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off-diagonal UGDFs related to both nucleons. The vertex factor V^^^li^lij I2) = 
^Xqxl{li^Q2,ki,k2) is the production amplitude of a pair of massive quark q and antiquark 
q with helicities Xq, Xq and momenta ki, k2, respectively. The longitudinal momentum 
fractions of active gluons are calculated based on kinematical variables of outgoing quark 
and antiquark: Xi = ^exp(+?/q) + ^exp{+yq) and X2 = ■^exp(-?/q) + ^exp(-?/q), 
where mg^t and rriq^t are transverse masses of the quark and antiquark, respectively, and i/q 
and Hg are corresponding rapidities. 

The off-diagonal UGDFs are written as 13 



fg {x',Xi^2,qi,2vqo,t,f^F) - Rg fgi^i,2,qi,2vf^F)^ (6-1) 

where Rg ~ 1.2 accounts for the single logQ^ skewed effect 3l[. The factor Rg here cannot 
be calculated from first principles in the most general case of off-diagonal UGDFs. It can 
be estimated only in the case of off-diagonal collinear PDFs when x' ^ x and xg = x~'^(l — 
x)" and then Rg = ^rix+^^ ■ considered kinematics the diagonal unintegrated 



densities can be written in terms of the conventional (integrated) densities xg{x, q^) as [30 



d I 

fg{x,qlfi^) = -^^[xg{x,q^t)^Tg{q^t.l^'^)] , (6.2) 

where Tg is the conventional Sudakov survival factor which suppresses real emissions from 
the active gluon during the evolution, so the rapidity gaps survive. 

In the framework of the /ct-factorization approach [32| the hard subprocess g*g* — ?■ qq 
gauge invariant amplitude y\^'x^{qi-iq2) reads 

V^lTMi, q2) = nlnZV^ll^^^\q,, g^, ^i, A;^), = (6.3) 
V^X''''{qi,q2) = -glY.i^'^^k\l)u,^{k,) X (t^jt-6^''(A;i,A;2) -t^,^/;^6'^'^(fc2,fci))t;A,-(fc2), 

i,k 

where Ep^cms = V^/^ is the c.m.s. proton energy, are the color group generators in 
the fundamental representation, u{ki) and f (/C2) are on-shell quark and antiquark spinors, 
respectively, b'^'^ and fe'^^ are the effective vertices arising from the Feynman rules in quasi- 
multi-Regge kinematics (QMRK) approach (33l |: 

^"''"'^-^ iq,-k,y-m^^ {k, + k2f ' ^^-^^ 
6-(A:„ k,) = 7-7^^^,^^,^7^ - 



{qi — ^2)^ — fn"^ {ki + h 



2 



2 



where T^^^{qi, ^2) is the effective three-gluon vertex. The effective ggg-vertices are canceled 
out when projecting the qq production amplitude Eq. fl6.3p onto the color singlet state. 
Since we will adopt the definition of gluon polarization vectors proportional to transverse 
momenta qi/2±, i-e. £1,2 ~ qi/2±/xi^2 (see below), then we must take into account the 
longitudinal momenta in the numerators of effective vertices (see Eq. 06.41) ). 

The SU(3) Clebsch-Gordan coefficient {3i,3k\l) = 5'^^ / y/Wl in Eq. (16. 3p projects out the 
color quantum numbers of the qq pair onto the color singlet state. Factor 1 / \/Wc provides 
the averaging of the matrix element squared over intermediate color states of quarks. 
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Therefore, we have the following amplitude 

(6.5) 

In the present calculations we use the GJR08 set of collinear gluon distributions [s^]. In 
the analogy to the CEP of Higgs boson, where renormalization and factorization scales are 
advocated to be /i^ = = = Mfj jssi], we apply the following prescription /i^ = M^g. 
Absorption corrections to the bare pp — ?• p{qq)p amplitude, which are necessary to be taken 
into account (to ensure exclusivity of the process), are included approximately by multiplying 
the cross section by the gap survival factors Sg = 0.1 for RHIC and Sq = 0.03 for the LHC 
energy. More details about exclusive production of heavy quarks can be found in our original 
paper [oj. Let us come now to presentation of our results. 

In FigJTS] we show rapidity distribution of c quarks from the exclusive mechanism shown 
in Fig. [T7] (solid line). We show the results for leading order (upper curves) and next-to- 
leading order collinear gluon distributions We observe large difference of results for LO 



and NLO gluon distribution especially at LHC. For comparison we show the contribution of 
central diffractive component discussed in section |Vl In this calculation we have included 
gap survival factors = 0.1 for y/s = 500 GeV and Sq = 0.03 for y/s = 14 TeV. The cross 
section for the exclusive mechanism is similar to that for the inclusive central diffractive 
mechanism. The exclusive production starts to dominate only at large c quark rapidities. 
Therefore a measurement of the cross section with double (both side) rapidity gaps may be 
not sufficient to single out the exclusive mechanism. Clearly other cuts would be necessary. 




FIG. 18: Distributions in rapidity of c quark/antiquark for the exclusive component at = 500 
GeV (left panel) and = 14 TeV (right panel). For comparison we show the central diffrac- 
tive contribution (dashed line). Different collinear gluon distributions were used to obtain the 
unintegrated gluon distribution according to the KMR prescription. 

Corresponding distributions in the c quark (c antiquark) transverse momentum are shown 
in FigJTOl The distribution for exclusive component extends to higher transverse momentum 
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than that for the central inclusive diffractive one. A lower cut on c quark (D meson) trans- 
verse momentum may therefore help to identify the exclusive component but will exclude 
a measurment of the integrated cross section for this component. A detailed Monte Carlo 
studies of final states of both components may help to find a better criterion to separate 
experimentally the two components. 




FIG. 19: Distributions in transverse momentum of c quark/antiquark for the exclusive component 
at ^/s = 500 GeV (left panel) and ^/s = 14 TeV (right panel) TeV. Different collinear gluon distri- 
butions were used to obtain the unintegrated gluon distribution according to the KMR prescription. 
For comparison we show the inclusive central diffractive contribution (dashed line). 



VII. CONCLUSIONS 

In the present paper we have calculated differential distributions for different subdominant 
contributions usually neglected in the literature when calculating production of cc pairs. 

Single and double photon induced processes are first class of mechanisms considered here. 
In calculating single particle distributions we have used a special set of parton distributions 
which includes photon as a parton of the proton. The calculation of the cross section is 
therefore very similar to that for the gluon-gluon fusion. The difference is only in color 
factors and a lack of the s-channel diagrams for photon induced processes. We have also 
included contributions when emitted photon which enters a hard process leaves a proton in a 
ground state. Those "elastic" mechanisms give similar contribution as the "deeply inelastic" 
mechanisms considered in the QCD-improved parton model. We have found that although 
individual contributions are very small, when added together they can give cross section of 
about 1 % of the inclusive one dominated by gluon-gluon fusion. In our analysis we have 
neglected resonance contributions when the photon leaves the remaining object in a proton 
excited state, e.g. in A(1220) resonance or other nucleon resonances. 

We have also discussed single and central diffractive production of cc pairs in the 
Ingelman-Schlein model. In these calculations we have included diffractive parton distri- 
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butions obtained by the HI collaboration at HERA and absorption effects neglected in some 
early calculations in the literature. The absorption effects which are responsible for the 
naive Regge factorization breaking cause that the cross section for diffractive processes is 
much smaller than that for the fully inclusive case, but could be measured at RHIC and 
LHC by imposing special condition on rapidity gaps. 

Finally we have discussed a fully exclusive diffractive production of cc. It was advocated 
recently that the cross section for this mechanism may be substantial. We have found here 
that both at RHIC and LHC its contribution is smaller than that for single diffractive one. 
In our opinion it is very timely to analyze if this contribution could be measured. This 
equires an analysis of the final state. We expect that the final state in single and exclusive 
production are different enough to set criteria to pin down the fully exclusive component. 
It is, however, not obvious if the central diffractive and purely exclusive mechanisms could 
be differentiated experimentally. They may look similar as far as rapidity gap structure is 
considered. We predict that the total contribution of central diffractive mechanism is similar 
to that for the exclusive one. In contrast the final state multiplicity can be expected to be 
different. A better analysis requires a Monte Carlo studies. 

We have not discussed an impact of diffractive mechanisms considered in the present 
paper on the fully inclusive cross section for cc pair production. This is a rather difficult 
task and goes beyond the scope of the present paper. 
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